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Abstract 

Manufacturing end-of-fiber optical components able to realize optical functions ranging from 
a simple lens to more complex functions such as mode selective components is a decisive but a 
priori complex task owing to the fiber core dimensions. Effective low cost methods allowing to 
grow polymer components by free-radical photopolymerization using the light coming out of the 
fiber have recently been reported. A novel phenomenological model of the photopolymerization 
process underlying is here given and used to simulate the polymer component growth in a three 
dimensional time-resolved manner. The simulations results are thus used to understand and 
optimize the component growth conditions, focusing particularly on the role of oxygen either 
present in the atmosphere or dissolved in the solution. 

PACS numbers: 060.2340 Fiber optics components; 160.5470 Polymers; 250.5460 Polymer waveguides-fibers 
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I. INTRODUCTION 



During the past decades, numerous methods for integrating microlenses on optical fibers 



(MOFs) have been reviewed. MOFs have been manufactured bv various methods such as 
ectric arc melting J| laser micromachining 2j, chemical etching 3] and photolithography 4, 



6[. Nevertheless, these procedures are time and energy consuming and imply stringent 



control of the experimental parameters. 



was recently reported 



8. 



A novel low cost manufacturing technique based on free-radical photopolymerization[7] 

10. 



ll|. It relies on photopolymerization within a photosensitive 



droplet attached to a cleaved fiber end thanks to surface tension. A beam of visible green 
light introduced at the other end of the fiber emerges beneath the droplet and thus initiates 



polymerization. 



As recent works show 12|, [l3|, [14j , light propagating in a photosensitive solution in which 



the refractive index of the polymer is higher than that of the initial formulation can ex- 



perience self-focusing and 



propagate in a spatial soliton manner as it is the case in other 



materials such as Kerr 15] or photorefractive media 







16l . Il7| . Therefore, since the polymer 



component remains firmly attached to the optical fiber thanks to hydrogen bonds, it can be 
inferred that the component grown following this procedure will be an extension of the fiber 
core and that its optical function will be governed by its index and shape. In the following, 
we will thus refer to it as polymer tip. 

In the following, a thorough 3D time resolved phenomenological model of this component 



growth on a cleaved fiber end is proposed. Starting with a review of the experimental 
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procedure, an abstract model able to describe the local polymerization rate will be presented, 
eventually leading to the growth modelling thanks to the coupling with a standard Beam 



Propagation Method [18l Il9l l20l| . It is thus followed by the analysis of the polymer tip growth 
conditions through simulation in order to understand and master the tip shape and optical 
function. The model presented in this paper is based on previously published preliminary 
work 8J, |9| and extends the previous models so as to allow their quantitative testing again 
experimental data. 

II. EXPERIMENTAL DESCRIPTION 
A. Photosensitivive solution 

1. Composition 

The photopolymerizable solution used in this work was introduced in previous papers 



21|,|22j]. It is made of three basic components: a sensitizer dye, an amine cosynergist, and a 
multifunctional acrylate monomer, pentaerythritoltriacrylate (PETIA) used as received from 
the supplier. After irradiation, this component forms the backbone of the polymer network. 
The cosynergist was methyldiethanolamine (MDEA), while eosin Y (2',4',5',7'-tetrabromo- 
fluorescein disodium salt) was used as sensitizer dye. This system was developed mainly 
because of its high sensitivity in the spectral region going from 450 to 550 nm (maximum at 
530 nm). In particular, this formulation is specially suited to be irradiated by an argon ion 
laser (514 or 488 nm), a frequency doubled Nd/YAG laser (532 nm) or a green He-Ne laser 
at 543.5 nm which we have used in the experiments described in this paper. In addition, 



this liquid system is very flexible, since it is possible to modify the nature or the proportions 
of the different components in order to adjust the physical and chemical properties of the 
formulation (viscosity, polymerization threshold energy, etc.). Results reported in this paper 
were obtained with mixtures containing 3 wt.% eosin and 8 wt.% methyldiethanolamine. 



2. Photopolymerization reaction 

After absorption of actinic light by eosin, the triplet state of the dye reacts with the 
amine to form radicals. Radicals initiate the polymerization of the monomer. Owing to the 
monomer multifunctionality, the polymer quickly develops into a 3D network. The nonpoly- 
merized formulation stays liquid and can be easily washed out with solvent, such as ethanol, 
leaving the polymerized part as a quasi- homogeneous transparent polymer component. 

However, as previous studies of this photopolymerization process have shown 23], oxygen 
dissolved within the formulation plays an inhibiting role in the photo-initiation process. As 
we will see in the following, the oxygen inhibition activity plays a determinant role in the 
polymer tip growth and thus largely governs its shape and optical function. 



B. Experimental set-up 

As is depicted in Fig[Q the experimental bench set up for the growth of the polymer tip is 
fairly simple and thus quite low cost. It basically consists in dipping a cleaved optical fiber 
in a photosensitive solution such as the one described above so as to obtain a hemispherical 
droplet attached to the end of the fiber by mere surface tension. Coupling a green laser light 
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(for instance a He-Ne laser at a wavelength of 543.5 nm) to the other end of the fiber allows 
light to emerge below the droplet and initiate polymerization. 

The rest of the process consists in adjusting laser power to the desired value and in letting 
the polymerization occur during the desired exposure time. The resulting polymer tip is 
still embedded in a non-polymerized monomer solution that can be cleansed by ethanol, for 
instance, thus leaving the polymer tip standing alone firmly attached to the fiber end. 

C. Sample results and applications 

An example of the obtained polymer tip is shown on Fig|2]with a detail shown on Fig|3J 
As expected, the tip looks like an extension of the fiber core and has a length roughly equal 
to the thickness of the surface tension driven droplet in which its growth took place. On 
the contrary, the curvature radius of its end is significantly smaller than that of the droplet. 
Indeed, if they were equal, the tip end would, at that scale, seem rigorously flat. This 
discrepancy will be thoroughly investigated in the following. 

Fig0 shows the polymer tip grown from a LP21 fiber mode. As can be seen, its shape is a 
real 3D mold of the electromagnetic intensity emerging from the optical fiber. Preliminary 
experiments show that this particular tip is a mode selective component: if a plane wave is 
shone onto it, the excited mode within the fiber happens to be the LP 2 i mode which created 
the tip. 

Other applications based either on the tip optical or mechanical properties can be fore- 
seen. As will be published elsewhere, the polymer tip is able to enhance laser diode to 
fiber coupling in a significant way. Furthermore, the polymer nature of this tip makes it an 
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attractive low cost non breakable tip for Scanning Near Field Optical Microscopy. 



III. PHENOMENOLOGICAL MODEL 

A. Photopolymerization rate empirical model 

The first step of the photopolymerization consists in a photosensitization of the dye - 
eosin in our case — which is excited in its triplet state by the incoming light. This reaction 
is followed by radical formation which consists in an electron transfer from the amine to the 
dye, followed by proton transfer. The polymerization is thus initiated and can propagate 
into a 3D polymer network which tends to solidify the solution. 

This reaction main pathway has however to compete with two side reactions involving 
O2 which thus inhibits the polymerization process. Indeed, O2 is able to deactivate the 
dye triplet state, thus impairing polymerization. It can also peroxydate the free radicals 

n 

and thus rapidly terminate the polymer chain[23j. Furthermore, these two reactions have 



been found |2^| to be much more probable than the polymerization process. Therefore, the 
deactivation reactions will dominate the polymerization process until O2 has disappeared 
from the solution. This is only a very crude model since it does not take into account O2 
diffusion within the solution towards O2 depleted region. However, this local model can be 
a good start because it can be expressed in terms of threshold energy. 

The threshold energy is in this case an energy density above which the time integrated 
optical energy density at a given point in space — which we will call the accumulated energy 
E acc in the following — has to grow before polymerization can start. It is directly linked 



to the concentration of dissolved O2 in the formulation to which it is roughly proportional. 
Indeed, a typical polymerization response curve is shown on FiglSJand exhibits the threshold 
energy E s . 



B. Complete 3D model of light propagation 

In view of this simple model, the simulation of the growth of the polymer tip requires on 
the one hand the measurement of the response curve as is given in FigEl and on the other 
hand the calculation of the accumulated energy E acc for every point in space within the 



25 
20 



and 



solution. The latter can be done knowing the electromagnetic field on the fiber enc 
using a standard Beam Propagation Method (BPM) based on Fourier Transforms jl8. 
The former allows to get access to the values of the threshold energy E s , the maximum 
index variation dn and to the polymerization rate linked to the derivative of the curve of 
FigElwith respect to the accumulated energy. Measurement methods can be found in the 



litterature 23|, |27( . As will be presented in section III, the phenomenological parameters not 
directly accessible via measurements (such as C s and C smax of equation |3] ) are going to be 
deduced by fitting the simulations to experimental results. 

The accumulated energy E acc is then computed iteratively using the BPM to know the 
energy density. The polymerization response curve is then used to evaluate the change in 
the refraction index at every point in space, thus providing data for the next BPM which 
will correspond to the next step in time. 

As the polymerization response curve used in the simulation is the results of experimental 
measurements, we have introduced an empirical model which allows its description as a 



function of the threshold energy E s , the polymerization rate immediately after the threshold 
— related to a — and the maximum index variation dn: 

An = dn x tanh ^ acc _ ^ s ^j x h (E acc — E s ) (1) 

where h is the Heaviside unit step function. 

Let us point out here that we have taken into account the fact that the polymerization 
efficiency does depend on the instant local intensity by adding a multiplier (3 to the local 
intensity whose value is taken between and 1 corresponding respectively to null efficiency 
and maximum one. However, at the typical intensities we use — 2,5mW within the fiber 



core — , the value of f3 is saturated to 1 23]. 

After a number of steps corresponding to the desired exposure time, this process yields 

a 3D matrix describing the refraction index variation in space. The shape of the polymer 

tip which is obtained after cleansing is then described by the polymer refraction index n p ( 
.51 < n p < 1.52 after irradiation, the index of the photopolymerizable being around 1.48 
2^). The shape of the tip is thus given by the surface corresponding to a refraction index 

constant value equal to n p . The typical tip shape obtained with this procedure is shown on 

Fig© 

The curvature of the hemispherical droplet is not taken into account in the simulation 
process. This is justified by the fact that the radius of curvature of the droplet (larger than 
the 125 microns fiber diameter) is much larger than the simulation window (on the order 
of fiber core, a few microns). Therefore, the hemispherical droplet is assumed flat in the 
simulation region. Furthermore, the droplet thickness of about 35 microns in the simulation 



region has been measured through electron beam microscopy. 

One particular aspect of the above described simulation scheme is that it is resolved in 
both 3D space and time. This means that light density and thus polymerization degree is al- 
lowed to be non-constant through time and non-uniform through space. Particularly, Figure 
6, as all simulation figures, has been computed while taking into account the inhomogeneous 
intensity distribution resulting from the propagation of a Gaussian beam. 

C. Dye photo-bleaching 

Our model also needs to consider the loss in absorption due to the consumption of absorb- 
ing dye that we must take into account in the BPM process. We have taken into account 
this phenomenon known as photobleaching in a similar way as the polymerization state. 
Experimental measurements allow to model phenomenologically the absorption represented 



unction of the accumulated energy E acc at every point 



by 7, the absorption coefficient, as a 
of the 3D space in the following way 



where d is the measurement sample thickness (i.e. an experimental constant), 70 the ini- 
tial absorption coefficient of the solution and k a characteristic energy density determined 
experimentally. 

Our simulations results show that this photobleaching curve has to be determined by a 
careful experiment as depicted in the litterature 



281 ] for it greatly influences the polymer tip 
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shape. In our case, the main parameter k has been measured to 12.5W/cm 2 . 
D. Oxygen diffusion within the solution 

A simple qualitative comparison between FigJ3] and FigO reveals a strong discrepancy 
between the experimental results and the simulation: the former show a rounded tip with a 
radius of curvature that goes down to the wavelength in particular experimental conditions 
whereas the latter evidences a flat end polymer tip. 

This discrepancy cannot be attributed to the curvature of the solution droplet because 
it is several orders of magnitude larger so that the polymer tip would indeed look flat if 
the curvature was only due to surface tension. Its origin has to be looked for into what the 
model does not account for: O2 diffusion into and within the droplet. 

The major source of O2 is the growth atmosphere which is generally the ambient atmo- 
sphere. Diffusion has then to be accounted for from the droplet surface into the solution for 
a few tenth of microns of length during the exposure time. The model presented here does 
not introduce the classical diffusion Fick law to model O2 diffusion for two reasons. The first 
one is a question of computational complexity, which would grow too high for a reasonable 
3D model to be developed. The second one is that the corresponding equation is too com- 
plicated. The spatially inhomogeneous incident intensity distribution in the formulation at 
the end of the fiber gives rise to a corresponding nonuniform distribution of both the viscos- 
ity and the reactives and photoproducts in the droplet. In particular, the inhomogeneities 
in the oxygen composition of the droplet should be taken into account, owing to different 
reactions with transient species and, simultaneously, to mass diffusion. 
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It is however possible to infer a phenomenological model from the experiments. Con- 
sidering that the main process is the diffusion of O2 from an infinite source into a medium 
with a homogeneous viscosity constant, the 2 concentration can be expected to follow an 
exponential law from a given concentration at the droplet surface. The consequence on 
the above described model is a spatial modulation of the threshold energy E s following 2 
concentration. This modulation can be described as follows: 



where C smax is linked to 2 atmospheric concentration on the droplet surface and D di ff is 
called the diffusion length, d is the droplet thickness and z the position of the current point 
in space counted from the fiber end. 

The same analysis can be carried out in the direction transverse to the beam propagation 
direction: when the beam propagates and initiates polymerization, it consumes oxygen. 
Following Fick's diffusion law, the 2 concentration gradient between the beam path and 
outside the beam path leads to an inward 2 flux towards the center of the beam. This will 
impair polymerization on the beam wings. We propose to model this the same way with a 
spatial modulation in the threshold energy as follows: 



where C smaX 2 is linked to 2 concentration within the droplet and -Dd«//2 is a diffusion length 
- a priori equal to D di ff. d 2 is the beam width on the order of the core radius and r the 
position of the current point in space counted from the beam axis. This phenomenological 




(3) 




(4) 
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model can be understood easily if the beam is considered to be homogeneous with a radius 
equal to di. However, if one needs to account for the Gaussian shape of the beam, this simple 
model needs to be adapted to account for the fact that 2 depletion varies continuously, 
potential leading to a raise in the threshold energy in the center of the beam by opposition 
to its wings. We will not account for that in this paper because we believe that the number 
of parameters arbitrarily introduced will grow too high for any phenomenological model 
to be viable. We will furthermore show that this refinement is not needed to explain the 
experimental measures. Consequently, undertaking this refinement would imply to reject 
the phenomenological approach and its computational simplicity and to solve directly the 
classical laws of diffusion. 

If we call C s the constant threshold energy due to dissolved 2 during the preparation 
process and taking into account the diffusion diffusion process modeled through expression 
EH and El the threshold energy can be spatially modeled as follows: 



where C smax and C smaX 2 are taken to have the dimensions of an energy density (or fluence). 
E. Sample results and model validation 

With this phenomenological model of the photopolymerization process, it it is possible to 
simulate the time-resolved growth of the polymer micro-tip. Figure U\ shows the computed 
growth of a polymer tip from the fundamental mode of a standard monomomode telecom 
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fiber excited at 543.5/im , to be compared to figure El which exhibits the experimental shape 
obtained. The beam profile shown on figure [7| does not seem to exhibit self-focusing as was 
discussed earlier — save for the last picture corresponding to a 5s exposure time. This is 
simply due to the diffraction length of the beam emerging from the fiber, which largely 
surpasses the droplet thickness — 117 jum vs. 35/zm . 

Figure |H1 shows the same build-up from a beam emerging from a monomode fiber at 
543.5/im . The diffraction length is in this case half the droplet thickness and the beam can 
easily be seen self-focusing, even exhibiting an oscillating behavior as predicted earlier 

IV. GROWTH CONDITIONS DESIGN 

A. Solution manufacturing and growth atmosphere control 

Figures El to [TT] alongside tables HI to IIHI show the simulations results for varying C s , C smax 
and Ddiff for a given solution and exposure time at a constant optical power. These figures 
allow to evidence the influences of these model parameters. One can read them easily once 
knowing that C s increases in the horizontal direction and C smax in the vertical one. 

The atmosphere in which the photosensitive solution is prepared associated to its sol- 
ubility in the solution plays a crucial role in the concentration of the oxygen dissolved 
homogeneously within the solution, the model parameter accounting for this phenomenon 
being C s . The analysis of the figures allows to show that an increase in the dissolved oxygen 
concentration leads to a decrease of the size of the tip. More precisely, it is the diameter of 
the polymer tip at the point where it is attached to the fiber which is mostly affected. We 
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can then infer that C s will play a role in the size of the modes excited in the tip. 

Furthermore, the influence of the oxygen diffusing from the atmosphere can also be re- 
trieved from these figures. An increase of C smax is shown to imply a decrease in the radius of 
curvature of the polymer tip end and thus a drastic change on the optical function provided 
by the tip. Let us point out here that we evidenced that it is the growth atmosphere oxygen 
concentration which designs here the polymer tip optical function. 

B. Solution viscosity and oxygen diffusion 

However, once again, the oxygen diffusion from the growth atmosphere towards the so- 
lution not only depends on the oxygen partial pressure of the atmosphere and the oxygen 
solubility but also on the solution itself and its ability to let the oxygen diffuse. Our model 
accounts for this phenomenon through the oxygen diffusion lengths D di fj and D di ff 2 . A 
careful analysis of figures El to [TT] and the corresponding tables allows to evidence that an 
increase of the oxygen diffusion length leads to a decrease in the radius of curvature of the 
tip end. 

Indeed, a rapid glance at the figures would make the reader think exactly the opposite 
but attention has to be put on the values of C s and C smax for each simulation. These 
varying concentrations of oxygen are simply due to the fact that, experimentally, the oxygen 
diffusing from the atmosphere does diffuse only in a thickness of a few microns and thus 
does not influence the photopolymerization process near to the fiber end. This leads to 
the following relationship which has to be verified for our model to be physically realistic: 



15 



C s C smax exp ( n ) . Furthermore, the homogeneously dissolved oxygen does not play 



any role at all near the interface with the oxygen rich growth atmosphere, which implies 
C S max ^> C s , thus leading to an overall relation which has to be verified for all our simulations 



A variation of Ddiff thus implies a correlated variation of C s and C smax . This, of course, 
does not imply physically that the oxygen diffusion length in the solution has an influence 
on the oxygen concentration in the atmosphere but rather that only a small range of varying 
Ddiff can physically correspond to a given oxygen concentration. Therefore, equationElgives 
us useful hints on which numerical values to use for our model parameters in order for it to 
describe correctly the experiments. 

The process of diffusion of O2 within the solution from dark regions to the bright ones 
which are 2 depleted is taken in account by the couple of parameters C smax2 and D di ff 2 . 
Their influence on the tip shape can be inferred and is confirmed by figurdT^l Raising either 
C sma x2 or D^iff 2 leads to a global narrowing of the tip by raising the threshold energy on 
its sides. 

C. Curvature radius prediction and model validation 

The model that is proposed in the previous sections is based on experimental data as 
far as both polymerization as a function of accumulated energy and the photo-bleaching of 




(6) 
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the dye are concerned. The four simulation parameters accounting for O2 diffusion into and 
within the formulation are on the contrary to be determined by iterative simulation and 
systematic comparison to experimental data. 

It is therefore a systematic experimental measurement campaign that has to be under- 
taken if one wishes to set the quantitative values for the model parameters. The experimental 
fit has then to be realized on a value that is both significant and quantitatively within reach 
both experimentally and through simulation. In this view we have chosen to measure in 
both cases the radius of curvature of the tip end. It is indeed the tip end which will largely 
govern the optical function achieved. 

In our case, the theoretical fit could not be automated because of the large amount of 
computing time it would have required: on the contrary, we have manually derived a set 
of parameters that allows a correct fit between theory and experiment, as can be judged 
by two significant experimental series which are reported on figures IT31 and ITU The former 
reports the tip end curvature radius as a function of the exposure time whereas the latter 
does as a function of the beam power that is injected in the fiber. 

Indeed, on the same figures are shown the simulation results obtained from the model 
parameters given in the caption. As can be expected, threshold energy variation due to the 
diffusion of the oxygen from the growth atmosphere — uncontrolled in the experiment — is 
much larger than that which is due to the diffusion of oxygen within the droplet. Further- 
more, the diffusion length found for both this phenomenon are equal. This again is expected 
since diffusion occurs in the same medium for both phenomenon. Moreover, the diffusion 
length found is around 0.1 /im . This value is coherent with previous results. Indeed, dif- 
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fusion coefficients of dye molecules in a polymer host were measured via the holographic 
grating relaxation techniques 29. I30I ]. A holographic grating was induced by crossing two 
equal intensity coherent beams issued from a pulsed laser in order to make the measurement 
of the slow translational diffusion coefficient feasible. After excitation, the dye concentration 
was gradually smeared out by diffusion. The corresponding optical grating relaxed with a 
time r = characteristic of the diffusion coefficient D, A being the fringe spacing. Taking 
10~ 7 cm 2 /s for the oxygen diffusion coefficient (by use of the Wilke-Chang correlation 3l| ) 
and A = 10/im (i.e. the lateral dimension of the irradiated zone corresponding to a cylin- 
drical symmetry), a diffusion time of 0.03 s is obtained, corresponding to a value lower than 
the exposure duration (a few seconds). Then, diffusion of oxygen from the non- irradiated 
part of the droplet towards the bright area where the monomer is converted into polymer 
is reasonable. A diffusion length of 0.1 /im is here self-consistent due to the reaction of 
oxygen with all the different active species present in the living polymer. The probability of 
reaction with an initiating radical, a living monomer radical or a monomer radical trapped 
in the crosslinked network is sufficiently high to obtain finally a diffusion length around 0.1 
fim during the growth of a polymer micro-tip. 

In this view, we believe that the quantitative model issued from this experimental fit is 
then able to predict the polymer tip shape for this photosensitive solution and for various 
experimental conditions. 

Therefore, this quantitative model can now be used to explain and interpret experimental 
data. For instance, both figures IT3~1 and ITU show a change of behavior after an input fluence 
of 10/xJ, the former showing a leap in the curve slope, the latter showing a large break in 
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the curve. The reasons for these behaviors is yet unclear but our simulations allow to show 
clearly that these behavior changes are related to the tip length. Indeed, the change occurs 
in both cases when the tip end reaches within an oxygen diffusion length from the droplet 
end: prior to the break, no influence from the growth atmosphere can be seen whereas after 
the break, it is atmospheric oxygen which rules the tip end shape. 

V. CONCLUSION 

The phenomenological model proposed in this paper is based on both experimental mea- 
surements and adaptive fit of the four parameters accounting for oxygen diffusion. The 
set of quantitative experimental measurements presented here allows to get the value for 
the model parameters for the photosensitive solution we used. Of course, using a different 
solution would imply another determination of the model parameters. However, once the 
parameters are set, our model allows to predict the polymer tip shape for various incident 
beam powers and exposure times. Once this base set, future work sould allow to estab- 
lish a direct link between our phenomenological parameters and the physical and chemical 
parameters characterizing the photopolymerization. 

We have not, at this point yet, predicted the optical function of the polymer tip. However, 
we can assume that the polymer tip, once cleansed with ethanol, shows an homogeneous 
refraction index equal to that of the polymer at its maximum of polymerization. Indeed, 
the cleansing will remove the polymer which has a refraction index lower than the polymer 
index n p = 1.51. The ambient light will then allow to reduce the slight index variations 
which may be left in the tip by simply allowing the completion of polymerization. 
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In conclusion, our model allows the prediction off the shape of an optical object which 
shows an homogeneous refraction index in which electromagnetic waves propagate linearly. 
We can thus expect to find some of the shelf software able to do the simulation. Figure 
ITo1 shows the results obtained for a simplified tip by the software FIMMPROP-3D from the 
society Photon Design which we thank by the way for the figure. 

The combination of our model and FIMMPROP-3D thus provides the tool to predict 
the tip optical function and thus allows the micro-optical manufacturer to grow end-of-fiber 
polymer micro-components at ultra low cost. 

The authors wish to thank the CNRS and Lorraine Region for their partial support of 
this work. 
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TABLE I: Simulation values for FigEI All simulations are done with an exposure time corre- 
sponding to Is. 
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TABLE II: Simulation values for Fig llOl All simulations are done with an exposure time corre- 
sponding to Is. 
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TABLE III: Simulation values for Fig llll All simulations are done with an exposure time 
sponding to Is. 
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List of figure captions 

1. Experimental setup: He-Ne 543. 5nm laser light emerges from the fiber into the surface 
tension driven photosensitive droplet, thus initiating the polymerization reaction. 

2. Electron micrograph of a 30 /im long polymer tip grown on the end of monomode fiber 
at 543. 5nm. Its basis width is on the order of the fiber core: 3 microns. 

3. Detail of the end of the polymer tip shown of FigJ21 

4. Polymer tip grown with a LP21 mode on the end of a standard telecommunication 
fiber. 

5. Typical polymerization curve expressed in terms of index of refraction variation as a 
function of the accumulated energy E acc , the maximum index variation between the 
monomer and the polymer being dn = 0.04. The existence of a threshold energy E s 
below which no polymerization can occur is evidenced. 

6. Typical shape of a 30/im long polymer tip computed using standard BPM and the 
accumulated energy model shown on Fig|S] without taking into account O2 diffusion 
into the solution. 

7. Simulated growth of a polymer tip on the cleaved end of monomode telecom 1.55/im 
fiber excited in 543.5/im in its fundamental mode (bottom), alongside the beam profile 
during the growth (top). The beam emerges from the fiber on the left. 

8. Simulated growth of a polymer tip on the cleaved end of monomode 543. 5nm fiber 
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(bottom), alongside the beam profile during the growth (top). The beam emerges 
from the fiber on the left. 

9. Simulated influence of dissolved oxygen (C a ) and atmosphere 2 concentration (C smax ) 
for a diffusion length D di jf = \\ivn. From (a) to (e) C s = 10 3 J/m 2 and C smax rises. 
From (f) to (j) C s = 10 4 J/m 2 and C smax rises. From (k) to (o) Cs = 10 5 J/m 2 and 
Cgmax rises. Precise values are given in Table UJ The fiber-end to which is attached 
the tip is located on the left-hand side. The hollow ends of tips that can be seen on 
figures 19*1 through ITT1 are a rendering artifact and mean that the tips are actually flat. 

10. Simulated influence of dissolved oxygen (C s ) and atmosphere O2 concentration {C smax ) 
for a diffusion length D di ff = 2/im. From (a) to (e) C s = 10 3 J/m 2 and C smax rises. 
From (f) to (j) C s = 10 4 J/m 2 and C smax rises. From (k) to (o) Cs = 10 5 J/m 2 and 
Cgmax rises. Precise values are given in Table HT1 The fiber-end to which is attached 
the tip is located on the left-hand side. 

11. Simulated influence of dissolved oxygen (C s ) and atmosphere O2 concentration (C sm ax) 
for a diffusion length D d iff = 3fim. From (a) to (e) C s = 10 3 J/m 2 and C sm ,ax rises. 
From (f) to (j) C s = 10 4 J/m 2 and C sma x rises. From (k) to (o) Cs = 10 5 J/m 2 and 
Cgmax rises. Precise values are given in Table 11111 The fiber-end to which is attached 
the tip is located on the left-hand side. 

12. Left: simulated polymer tip without taking into account diffusion within the droplet. 
Right: droplet inner diffusion taken in account. Simulation data: C s =5x 10 2 J/m 2 , 
C smax =10 8 J/m 2 , D di ff 2 =D dif f =0.1/mi, C smax2 =10 4 J/m 2 and an exposure time of 
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Is. The dotted line around the right-hand tip is a shadow of the left-hand one meant 
to ease comparison. 

13. Experimentally measured curvature radii are compared to simulated ones for varying 
exposure times and for a beam power of 2,5/xW, an experimental eosin concentra- 
tion of 3% and the following simulation data: C s = 5xl0 2 J/m 2 ,C smax =10 8 J/m 2 , 
D di ff=D di ff 2 = 0,1/xm, C sma:K 2=10 4 J/m 2 and a = 10. 

14. Experimentally measured curvature radii are compared to simulated ones for varying 
output beam powers and for an exposure time of Is, an experimental eosin concen- 
tration of 3% and the following simulation data: C s = 5xl0 2 J/m 2 ,C sma x =10 8 J/m 2 , 
Drfj//=D rf jjj 2 = 0,1/iin, C smaa; 2=10 4 J/m 2 and a = 10. The straight lines are a mere 
guide to the eyes. 

15. Simulated propagation of light in a sample polymer tip: the initially plane wave un- 
dergoes strong focusing prior to reaching the tip end, thus providing high enlargement 
of the numerical aperture of the system. 
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FIG. 1: Experimental setup: He-Ne 543. 5nm laser light emerges from the fiber into the surface 
tension driven photosensitive droplet, thus initiating the polymerization reaction. 
2085 
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FIG. 2: Electron micrograph of a 30 /im long polymer tip grown on the end of monomode fiber 
at 543. 5nm. Its basis width is on the order of the fiber core: 3 microns. 
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FIG. 3: Detail of the end of the polymer tip shown of Fig[2l 
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FIG. 5: Typical polymerization curve expressed in terms of index of refraction variation as a 
function of the accumulated energy E acc , the maximum index variation between the monomer and 
the polymer being dn = 0.04. The existence of a threshold energy E s below which no polymerization 
can occur is evidenced. 
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FIG. 6: Typical shape of a 30/xm long polymer tip computed using standard BPM and the accu- 
mulated energy model shown on Fig|S] without taking into account O2 diffusion into the solution. 
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Exposure time : 5s 




FIG. 7: Simulated growth of a polymer tip on the cleaved end of monomode telecom 1.55/xm 
fiber excited in 543.5/mi in its fundamental mode (bottom), alongside the beam profile during the 
growth (top). The beam emerges from the fiber on the left. 
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FIG. 8: Simulated growth of a polymer tip on the cleaved end of monomode 543. 5nm fiber 
(bottom), alongside the beam profile during the growth (top). The beam emerges from the fiber 
on the left. 
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FIG. 9: Simulated influence of dissolved oxygen (C s ) and atmosphere O2 concentration (C smax ) 
for a diffusion length Ddiff = 1/um. From (a) to (e) C s = 10 3 J/m 2 and C smax rises. From (f) to 
(j) C s = W 4 J/m 2 and C smax rises. From (k) to (o) Cs = W 5 J/m 2 and C smax rises. Precise values 
are given in Table HJ The fiber-end to which is attached the tip is located on the left-hand side. 
The hollow ends of tips that can be seen on figures IU1 through II II are a rendering artifact and mean 
that the tips are actually flat. 
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FIG. 10: Simulated influence of dissolved oxygen (C s ) and atmosphere O2 concentration (C smax ) 
for a diffusion length D^iff = 2[/,m. From (a) to (e) C s = 10 3 J/m 2 and C smax rises. From (f) to 
(j) C s = W A J/m 2 and C smax rises. From (k) to (o) Cs = W 5 J/m 2 and C smax rises. Precise values 
are given in Table ITII The fiber-end to which is attached the tip is located on the left-hand side. 
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FIG. 11: Simulated influence of dissolved oxygen (C s ) and atmosphere O2 concentration (C smax ) 
for a diffusion length Ddiff = 3^m. From (a) to (e) C s = 10 3 J/m 2 and C smax rises. From (f) to 
(j) C s = W 4 J/m 2 and C smax rises. From (k) to (o) Cs = W 5 J/m 2 and C smax rises. Precise values 
are given in Table ITTT1 The fiber-end to which is attached the tip is located on the left-hand side. 
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FIG. 12: Left: simulated polymer tip without taking into account diffusion within the 
droplet. Right: droplet inner diffusion taken in account. Simulation data: C s =5xl0 2 J/m 2 , 
C sm ax=W 8 J/m 2 , D di ff 2 =D di ff =0.1/xm, C smaX 2=W 4 J /m 2 and an exposure time of Is. The 
dotted line around the right-hand tip is a shadow of the left-hand one meant to ease comparison. 
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FIG. 13: Experimentally measured curvature radii are compared to simulated ones for vary- 
ing exposure times and for a beam power of 2,5/iW, an experimental eosin concentration of 3% 
and the following simulation data: C s = 5xl0 2 J/m 2 ,C smax =10 8 J/m 2 , T>aiff=^diff2= 0,1/zm, 
C sma x2=10 4 J/m 2 and a = 10. 
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FIG. 14: Experimentally measured curvature radii are compared to simulated ones for varying 
output beam powers and for an exposure time of Is, an experimental eosin concentration of 3% 
and the following simulation data: C s = 5xl0 2 J/m 2 ,C smax =10 8 J/m 2 , T>diff=^diff2= 0,1/xm, 



C sma x2=10 J/m 2 and a = 10. The straight lines are a mere guide to the eyes. 
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FIG. 15: Simulated propagation of light in a sample polymer tip: the initially plane wave undergoes 
strong focusing prior to reaching the tip end, thus providing high enlargement of the numerical 
aperture of the system. 
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